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ABSTRACT 

Cooling flows in galaxy clusters and isolated elliptical galaxies are a source of mass for fueling accretion 
onto a central supermassive black hole. We calculate the dynamics of accreting matter in the combined 
gravitational potential of a host galaxy and a central black hole assuming a steady state, spherically 
symmetric flow (i.e., no angular momentum). The global dynamics depends primarily on the accretion 
rate. For large accretion rates, no simple, smooth transition between a cooling flow and an accretion flow 
is possible; the gas cools towards zero temperature just inside its sonic radius, which lies well outside the 
region where the gravitational influence of the central black hole is important. For accretion rates below 
a critical value, however, the accreting gas evolves smoothly from a radiatively driven cooling flow at 
large radii to a nearly adiabatic (Bondi) flow at small radii. We argue that this is the relevant parameter 
regime for most observed cooling flows. The transition from the cooling flow to the accretion flow should 
be observable in M87 with the Chandra X-ray Observatory. 



Subject Headings: galaxies: cooling flows - accretion, accretion disks - galaxies: individual (M87) 



1. INTRODUCTION 

X-ray observations of clusters of galaxies and isolated 
elliptical galaxies indicate the presence of significant quan- 
tities of hot gas whose radiative cooling time is much less 
than a Hubble time (see, e.g., Fabian 1994 for a review). 
This is generally interpreted as evidence for a cooling flow, 
i.e., the radiatively driven accretion of gas onto a central 
cluster galaxy or towards the center of the host galaxy 
(e.g., Cowie & Binney 1977; Fabian & Nulsen 1977; Math- 
ews & Bregman 1978). The inferred accretion rates range 
from ~ O.IM0 yr~^ to ~ lO^M© yr~^ and are often in- 
ferred to decrease substantially with decreasing radius, 
roughly as M oc r (e.g., Stewart et al. 1984). The decreas- 
ing M is probably evidence for an inhomogeneous flow, in 
which some of the gas (that which is denser) cools out of 
the hot flow at large radii. 

It has long been recognized that a cooling flow may be 
an important source of fuel for accretion onto a central su- 
permassive black hole in the host galaxy (e.g.. Bailey 1980, 
Nulsen, Stewart, & Fabian 1984). The goal of this paper is 
to calculate in some detail the dynamics of accreting mat- 
ter in the central regions of a galaxy which harbors both 
a cooling flow and a supermassive black hole. We restrict 
ourselves to spherically symmetric flows. 

Our model and its assumptions are discussed in the next 
section (§2). In §3 we present solutions to our model prob- 
lem and explain the physics which governs the dynamical 
transition from a cooling flow to an accretion flow. Finally, 
§4 sums up and discusses the observational implications of 
our work. 



2. THE MODEL 

We consider a simplified model which we believe cap- 
tures much - although not all (§4) - of the interesting 
physics of the transition from a cooling flow to an accre- 
tion flow. The primary assumptions of our model are: (1) 



steady state and (2) spherical symmetry. Aside from adia- 
batic compression, we assume negligible heating of the gas, 
and we approximate its cooling as due to bremsstrahlung 
emission. 

Cooling flows are inferred not to conserve mass. Theo- 
retically, this can be understood by noting that the cooling 
time of the gas is comparable to the inflow time. Over- 
dense "blobs" can therefore condense out of the flow be- 
fore accreting. We model this by including a mass drop 
out term in the continuity equation (cf White & Sarazin 
1987): 
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The quantities M, p, v, 0, c^, T, s, and 7 represent the 
mass accretion/inflow rate, the mass density, the radial ve- 
locity (negative inwards), the gravitational potential, the 
isothermal sound speed, the temperature, the entropy per 
unit mass, and the gas adiabatic index, respectively. The 
quantity B k, 5.6 x 10"*^^ ergs cm^ g^^ is the normalization 
for the bremsstrahlung emission. 

The efficiency of mass drop out depends on the ratio 
of the infiow time of the gas (tjnf = 'r/\v\) to the cooling 
time (tcooi = Cs/Bp). When t-„^f ^ tcooi (as in a cooling 
fiow), mass drop out can be important since there is time 
for overdense regions to condense out. When tint ^ tcooi, 
however (as in Bondi's spherical accretion solution), mass 
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drop out is negligible since the eas accretes before cooling. 
The parameter q in equation ffl) normalizes the efficiency 
of mass drop out; we expect g ~ 1. 

We take the gravitational potential to be given by 
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The left most term in equations and (||) is the grav- 
itational potential due to a black hole of mass M and 
Schwarzschild radius Vg = 2GAd/c^\ the l/(r — rg) mim- 
ics the effects of General Relativity (Paczyhski & Wiita 
1980). 

HST observations of the centers of elliptical galaxies in- 
dicate that those which harbor massive black holes have 
central surface brightness "cores." The surface brightness 
rises steeply with decreasing radius at large radii, but flat- 
tens out in the inner portions of the galaxy (e.g., Lauer 
et al. 1995; Faber et al. 1995; Kormendy & Richstone 
1995)_F] Our mass model for the galaxy in equations (jij) 
and is intended to reflect these observations. Outside 
a break radius r?,, the galaxy has a constant velocity dis- 
persion (cr) while inside that radius we assume that light 
traces mass (aside from the black hole) and so the enclosed 
mass profile flattens (/3 < 1) in accord with the flattening 
surface brightness. 

Typical observed values of (3 and are w 0.25 and w 1 
kpc, respectively (Faber et al. 1995); we use such values 
in our numerical results of §3.1. For analytical estimates, 
however, we take r?, — * 0, i.e., we model the galaxy as 
having a constant velocity dispersion everywhere. 

The gravitational potential is predominantly that of the 



galaxy for 
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black hole for r rtr, where the transition radius is (for 
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For a galaxy with cr w 300 km in the Virgo cluster, 
a distance « 20 Mpc away, rtr corresponds to « 0.5 arc- 
sec. This is comparable to the angular resolution of the 
Chandra X-ray Observatory (CXO), indicating that the 
presence of a central black hole in nearby ellipticals may 
have observable effects on the cooling fiow X-ray emission 
(see §4.1). 

3. MODEL SOLUTIONS 

We are interested in solutions to our model problem 
which undergo a subsonic to supersonic transition at a 
sonic radius r^. Rewriting equations (|l|)-(§) yields 
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and 
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Since Z? = at the sonic point, we must have iV = for a 
smooth transition. 

We take our outer boundary conditions to be specified 
values of po and Co, the density and sound speed, at an 
outer radius ro. Equations (|l|)-(^) have two eigenvalues, 
taken here to be the sonic radius, rg, and the accretion 
rate at the sonic radius, M{rs). We find our solutions by 
shooting out from r^ and adjusting M (r^) and rg to satisfy 
the outer boundary conditions. 

For most of this paper we scale our models to obser- 
vations of M87: ro « 100 kpc, po « lO'^"^ g cm-^, 
cr « 300 km s^^, Co ~ cr (e.g., Stewart et al. 1984) and 
Af « 3 X IO^Mq (Harms et al. 1994; Ford et al. 1994; 
Macchetto et al. 1997). In M87, M is inferred to decrease 
from « 10 Mq yr"i at w 70 kpc to ^ 1 Mq yr"^ at a few 
kilo-parsecs (Stewart et al. 1984). As shown below, this 
is reasonably well captured by a g = 0.6 model. 

Although strictly an eigenvalue of our problem, an ex- 
cellent estimate of the value of M at the outer bound- 
ary, M(ro) = Mo, can be made by imposing the re- 
quirement that the gas be a cooling flow at ro. This 
requires that the inflow time of the gas be comparable 
to the local radiative cooling time which in turn requires 
u(ro) = Vo ~ —Bporo/co and therefore (cf Fabian & 
Nulsen 1977) 
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The primary question is what happens to matter accreting 
at this rate at smaller radii. 

In the absence of a central point mass a cooling flow 
with Mo given by equation (|o|) would have the following 
structure: 

Cs~a, Af ocr39/(9+2), p oc r-3/(«+2) , & oc r(«-i)/(«+2) . 

This is readily verifled by substitution into equations (|1|)- 
(|^) under the assumption that the gas is subsonic and that 
the potential is logarithmic (M = 0, rf, — > 0). Note, how- 
ever, that if g < 1 the above scalings show that the Mach 
number of the flow increases as the radius decreases. At 
ro, the Mach number is 
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iNote, however, that the surface brightness does not become constant even on the smallest scales accessible by HST. It still increases with 
decreasing radius, but with a noticeably smaller power law slope. 
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Therefore, there is a critical radius at which the flow would 
make a sonic transition: 



r = r„it«roK/co)(«+2)/(i-«). 



(13) 



The key physics which determines the nature of the cool- 
ing flow to accretion flow transition is whether rtr ^ Tcrit or 
^ ^crit- That is, is the radius where the gravitational 
influence of the black hole becomes important inside or 
outside of the radius where, in the absence of the point 
mass, the cooling flow would have its sonic point? 

If Ttr ^ Tcrit the cooling flow is oblivious to the existence 
of the point mass. It undergoes its sonic transition with 
little modification because of the black hole. Inside of the 
sonic point, the gas quickly cools towards zero tempera- 
ture. The reason is that in the supersonic zone a cooling 
flow in a logarithmic potential has v oc [ln(rs/r)]"'^/^ and 
p (X r^'^v^^ ~ r~^. The velocity (density) therefore rises 
more slowly (quickly) with decreasing radius than in the 
subsonic zone, where the cooling time was of order the in- 
flow time. Consequently, in the supersonic zone the cool- 
ing time becomes much shorter than the inflow time and 
the gas cools rapidly (see, e.g.. Fig. 4 of Sarazin & White 
1987). In this case, a simple cooling flow to accretion flow 
transition is not possible. Whether such a cooling flow is 
an important source of mass for an accretion flow depends 
critically on the fate of the rapidly cooling matter inside 
the sonic point (e.g., does it form stars which lose mass to 
feed the central black hole?). 

For M87, rtr ~ 0.15 kpc; this is comparable to Tcrit if 
(7 = 0. Many other observed cooling flows have Mo 3> 
10 Mq yr~-^ (M87's value) and would have rcrit ^ rti if 
q were equal to 0. Significant mass drop out is, however, 
generally inferred in such cooling flows (for M87, q « 0.6 is 
appropriate; see below). From equation ( pi] ) it follows that 
if q 1 the Mach number of the cooling flow is roughly 
constant. This insures that the gas remains subsonic down 
to the radius where the potential of the black hole begins 



l/(r — Tg) potential]^ We have checked this carefully by 
numerically searching the Mo and space; we find that 
there is no transonic cooling flow solution for r < rtr. 

The upshot of this analysis is that the cooling flow makes 
a smooth transition to a nearly adiabatic Bondi flow in 
the vicinity of rtr- This provides a direct dynamical link 
between the cooling flow and the accretion flow onto the 
central black hole. 

The details of the flow structure near rtr are shown in 
Figure 1, where we plot the temperature, density, Mach 
number, and the ratio of the cooling time to the inflow 
time for models with q = 0.2, 0.6, and 1. Figure 2 shows 
the accretion rate as a function of radius for these mod- 
els. The boundary conditions at Tq = 100 kpc are those 
appropriate for M87. In M87, M is inferred to decrease 
from « 10 Mq yr"^ at « 70 kpc to 1 Mq yr"^ at a 
few kilo-parsecs (Stewart et al. 1984). This is reasonably 
well captured by our q = 0.6 model. The additional values 
of q are shown to indicate the range of expected behav- 
ior. Note that the models in Figures 1 and 2 assume that 
rh — 0, i.e., that the galaxy's potential is logarithmic ev- 
erywhere. Modifications due to alternative potentials are 
discussed in §3.1. 

The conceptually most important result in Figure 1 is 
the ratio of the cooling time to the infiow time. This is 
~ 1 at large radii where the accreting matter is indeed a 
cooling fiow but it increases rapidly at smaller radii as the 
fiow undergoes the transition to nearly adiabatic (Bondi) 
accretion. Associated with this transition is a decrease in 
the importance of mass drop out. In the Bondi regime, 
^cooi ^ tini and so overdense regions do not have time 
to condense out of the flow. This accounts for the nearly 
constant accretion rate for r ^ 300 pc in Figure 2. 



3.1. The Importance of the Potential 

To illustrate the impact of the underlying potential of 
the host galaxy on the structure of the cooling fiow, con- 
sider the analytical solution to equations (|])-(^ in the 
region where the galaxy's mass profile fiattens (i.e., for 
rtr ^ r < rt and (3 < 1): oc r'^+i, p oc r-^/'^-f/'^^ 
and V oc r^'^/'*"'"^/'* (for simplicity, we have taken q = 0). 
Because the cooling fiow is roughly virial, the gas temper- 
ature decreases inward. The Mach number of the resulting 
solution therefore increases rapidly inwards, as r'^/5/4-5/4_ 
For observed values of /3 « 0.25, the Mach number in- 
creases roughly as r~^, much more rapidly than the r~^/^ 
scaling in the logarithmic potential. One might worry that 
this would change the results of the previous section by 
forcing the cooling fiow to have its sonic point outside of 
rtr- Our calculations show, however, that this does not 
happen. 

Figure 3 shows the temperature, density, Mach number 
and ratio of cooling time to inflow time for three models of 
M87 with q — 0.6, and Figure 4 shows the corresponding 
predicted X-ray surface brightness proflles. The models in 
Figures 3 and 4 assume three different forms for the po- 
tential of the host galaxy: one purely logarithmic (solid 
line; the same model is also shown by a solid line in Figs. 
1 & 2) and two with flattening mass profiles; /? « 0.25 in- 

^ Since the Mach number of the cooling flow solution in the 1/r potential decreases with decreasing radius the only way it could ever make 
a sonic transition is if it were forced to by a singularity in the potential (or by some other boundary condition at small radii) . 



to dominate the potential of the galaxy (i.e., rtr <; '"crit)- 
We now turn our attention to this regime. 

There is one conceptual complication in understanding 
the cooling flow to accretion flow structure in the regime 
where rtr ^5 ^crit- Namely, in the point mass potential 
(i.e., for r ^ rtr), there are two viable analytical solutions 
to equations (0)-(||). One is Bondi's adiabatic spherical 
accretion solution (Bondi 1952). The other (valid only for 
7 > 11/7) is a cooling flow solution (cf Fabian & Nulsen 
1977) which has oc r~^, p oc r"^/* and v oc r^^^^. It is 
not a priori obvious which of these solutions the cooling 
fiow at large radii (in the galaxy's logarithmic potential) 
will match onto at small radii. One might in fact suspect, 
on the basis of "continuity" arguments, that the cooling 
fiow at large radii would match onto a cooling flow at small 
radii. 

Our calculations indicate, however, that this is not so. 
The reason is simply that there is no transonic cooling 
fiow solution in the r^^ potential. The Mach number 
of the cooling fiow solution in a point mass potential de- 
creases with decreasing radius as r^^'^; this solution cannot 
be matched to a sonic transition even with the introduc- 
tion of a singularity into the point mass problem via the 
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side Tjh « 0.7 kpc (dotted line) or r;, « 3 kpc (dashed-dot 
line)! 

As illustrated by Figures 3 and 4, the quantitative de- 
tails of the flow structure depend (not surprisingly) on 
the form of the potential on kilo-parsec scales. The most 
important result of this subsection, however, is that even 
when the mass profile flattens, the cooling flow still has 
its sonic point well inside r^^ and therefore still makes 
a smooth transition to a nearly adiabatic accretion flow. 
The qualitative picture outlined in the previous section is 
therefore unchanged. This is partially because the inferred 
Mach number in M87 at a few kilo-parsec is sufficiently 
small (w 0.04; see Fig. Ic) that even if the Mach num- 
ber increases linearly at small radii the black hole still be- 
comes important before the cooling flow has its sonic point. 
In addition, because the temperature of the gas decreases 
somewhat if the galaxy's potential flattens, the black hole 
becomes important for the dynamics of the gas at larger 
radii (the transition radius, rtr, is effectively oc cj^ when 
Cs < a). Note, for example, that fcooi/^inf is actually larger 
at smaller radii in the model with rf, = 3 kpc. This is 
because the decreasing temperature due to the flat mass 
profile implies that the flow undergoes its transition from 
a cooling flow to an accretion flow at a larger radius. 

4. DISCUSSION 

The theoretical picture which emerges from the above 
analysis is a very intuitive one: cooling flows in ellipti- 
cal galaxies remain subsonic down to roughly kilo-parsec 
scales at which point they are gravitationally "captured" 
by the central supermassive black hole. They then un- 
dergo a transition to a nearly adiabatic Bondi flow and 
accrete onto the black hole. There is consequently a direct 
and simple dynamical transition from a cooling flow to an 
accretion flow. 

Our calculations indicate that only if the cooling flow 
accretion rate on scales of a few kilo-parsecs is ^ 10 — 
100 Mq yr~^ will the cooling flow have its sonic transition 
before being captured by the black hole (the precise M de- 
pends on the black hole mass - assumed ^ lO^M© - and 
the galaxy potential, i.e., (3, a, and r^; less massive black 
holes, larger cr, and flatter galaxy mass profiles - larger 
and smaller P - require smaller M). In this case, the gas 
cools rapidly just inside its sonic point, which lies outside 
the region where the gravitational infiuence of the black 
hole is important. Although M <; 10 — 100 Mq yr^^ is 
readily achieved on ~ 100 kpc scales, mass drop out due 
to thermal instability, and the consequent decrease of M 
with decreasing radius, implies that observed cooling flows 
typically lie in the parameter regime where a smooth tran- 
sition from a cooling flow to a nearly adiabatic accretion 
flow is the viable transonic solution. We find that this con- 
clusion is valid even given uncertainties in the mass profile 
of the host galaxy on kilo-parsec scales (§3.1). 

It is interesting to note that all of the solutions shown 
in Figures 1-5 have Ts <^ rn-] that is, the sonic point lies 
well inside the radius where the potential changes from 
that of the galaxy to that of the black hole. Nonetheless, 
the dynamics of the accreting gas changes from that of a 
cooling flow to that of a Bondi flow at, or even outside of, 

^M87's surface brightness profile is observed to flatten around r;, « 



rtr . The requirement that the flow undergo a smooth sonic 
transition influences its structure out to radii <; rtr 3> r^. 

The primary assumption of our analysis is that there 
is negligible angular momentum in the inflowing material. 
Angular momentum may, however, become important on 
the small scales of interest here; if so, a transition to a 
Bondi flow will not be viable. In this case, we believe that 
the cooling flow will generally undergo a transition to an 
advection-dominated accretion flow (ADAF), the hot ac- 
cretion flow analogue of Bondi's solution for flows with an- 
gular momentum (Narayan & Yi 1994, 1995; Abramowicz 
et al. 1995). It is well known that ADAFs can only exist 
below a critical accretion rate, Mc- In the present context, 
this can be estimated by requiring the cooling time at rtr 
to be greater than the inflow time, which yields 

Me 10 Mq yr-1 ( -4^) ( ^ t) ■ (14) 

" VlO^o/ ^300 km s-i/ ^ ' 

In equation (|l^), a is the ratio of the inflow velocity to 
the local Keplerian speed, and is roughly the Shakura & 
Sunyaev (1973) viscosity parameter. 

The properties of the transition from a cooling flow to 
a Bondi flow described in §3 are determined primarily by 
the low radiative efficiency of the Bondi solution. Con- 
sequently, the qualitative picture outlined in this paper 
should apply equally well to the transition from a cooling 
flow to an ADAF, since the latter is also radiatively ineffi- 
cient (the quantitative details may, of course, change). If, 
however, the accretion rate supplied by the cooling flow 
to rtr is Mc, the cooling flow may not be able to evolve 
into an ADAF; it will instead likely collapse to a thin disk. 
We plan to explore this question in a subsequent paper. 

4.1. Observational Prospects 

To conclude, we discuss the observational prospects for 
observing the transition from a cooling flow to an accre- 
tion flow with the excellent (< 1 arcsec) angular resolution 
of the CXO. We focus our discussion on M87 in the Virgo 
cluster, as it is the closest large cooling flow. For a distance 
to M87 of « 18 Mpc, one arcsec is « 0.1 kpc. The basic 
inference from Figures 1 & 3 is that on scales of 0.1 kpc, 
the dynamics of accreting gas in M87 is far from that of 
a cooling flow. The ratio of the cooling time to the inflow 
time is likely « 10 or larger. 

One prediction of our model is that the inferred gas tem- 
perature in M87 should increase with decreasing radius on 
sufficiently small scales (when the black hole's potential 
dominates). In particular, all of our models show an in- 
creasing temperature profile for r <, 300 pc. Note however, 
that because the gas temperature in nearly all models is 
roughly virial (and therefore traces the potential), such 
a temperature increase would not necessarily confirm the 
cooling flow to accretion flow transition. It would merely 
show that the gas dynamics is dominated by a centrally 
condensed mass distribution (namely, the black hole). 

The second prediction of our model is that the X-ray 
surface brightness profile, while centrally peaked, should 
be suppressed with respect to that expected from a cool- 
ing fiow without a central point mass. This is because the 

7 kpc (Faber et al. 1995). 
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transition from a cooling flow to an accretion flow entails 
a decrease in radiative efiicicncy. A suppressed X-ray s\ir- 
face brightness profile is, however, also the characteristic 
signature of a cooling flow with mass drop out. 

At a minimum, our analysis indicates that signiflcant 
care must be taken in interpreting the X-ray surface 
brightness profiles of elliptical galaxies on 1 kpc scales. 
Standard analysis techniques, based on the assumption of 
a cooling flow, could incorrectly infer a radially decreas- 
ing accretion rate when in fact the accretion is via a Bondi 
flow with constant M and it is the radiative eSiciency that 
decreases. As shown in §3 (see Figs. 2 & 5), M is expected 
to be nearly constant on the small scales of interest here 
since the transition to a Bondi-like flow strongly suppresses 
mass drop out due to thermal instability. 

Given an observed X-ray luminosity Lx and sound 
speed Cg, the true local accretion of the flow is 

C2 iinf ^ ' 

Our calculations indicate that tcooi/^inf is likely to be ^ 10 
on « 0.1 kpc scales in M87 so that accretion rates inferred 
assuming a cooling flow (fcooi ~ Unf, M w Lx/c1) can be 



underestimates by up to an order of magnitude. This is 
shown explicitly in Figure 5 where we compare the true ac- 
cretion rate for our baseline (g = 0.6) model of M87 (solid 
line) with that which we would infer if we took the pre- 
dicted X-ray surface brightness and temperature profiles 
and assumed a pure cooling flow (dotted line) . 

X-ray spectra may break the degeneracy between a 
change in radiative efliciency and mass drop out. The pres- 
ence of co-spatial multi-temperature components in cool- 
ings flows is indicated by a detailed examination of their X- 
ray spectral lines (Canizares et al. 1979, 1982; Mushotzky 
et al. 1981). This is interpreted as evidence for mass 
drop out in the cooling flow. Consequently, the absence of 
evidence for such multi-temperature components (or a de- 
crease in their prominence) , together with a suppression in 
the X-ray surface brightness at small radii, would provide 
considerable support for our model of the cooling flow to 
accretion flow transition. 

We acknowledge support from an NSF Graduate Re- 
search Fellowship (EQ) and NSF Grant 9820686 (RN). We 
thank the referee for useful comments which signiflcantly 
improved this paper. 
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Fig. 1. — Numerical models of the transition from a cooling flow to an accretion flow. The temperature, density, ratio of 
cooling time to inflow time, and Mach number are shown (clockwise, from top left). The outer boundary conditions are 
chosen to represent M87. The q = 0.6 model reproduces M87's radially varying M reasonably well (Fig. 2). Models for 
two additional values of q are shown to indicate the range of expected behavior. 
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Fig. 2. — The accretion rate as a function of radius for the three models whose structure is shown in Figure 1. The 

q = 0.6 model roughly represents M87, for which observations show a decrease in M from « 10 yr"-'^ at « 70 kpc to 
< 1 Mq yr"-'^ at a few kilo-parsccs. The presence of a central black hole (taken to have M = 3 x 10^) and the associated 
transition to a nearly adiabatic accretion flow leads to a complete suppression of mass drop out on small scales. This 
accounts for the constant accretion rate on scales less than « 300 parsecs. 
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Fig. 3. — Numerical models of the cooling flow to accretion flow transition for three possible underlying potentials for the 
host galaxy (see eqs. Q and |^). The outer boundary conditions arc chosen to represent M87. The temperature, density, 
ratio of cooling time to inflow time, and Mach number are shown (clockwise, from top left). 
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4. — Predicted X-ray surface brightness profiles for the three models whose structure is shown in Figure 3. 
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Fig. 5. The solid lino shows the true radially varying accretion rate for our q = 0.6 model of M87 (pure logarithmic 
potential; solid line in Figs. 3 & 4). The dotted line shows the accretion rate inferred from the predicted X-ray surface 
brightness and temperature profiles if one assumes that the gas is a pure cooling flow. The true accretion rate is nearly 
constant on the small scales shown in this figure because the transition from a cooling flow to a nearly adiabatic accretion 
flow suppresses mass drop out due to thermal instability. The decreasing radiative eflaciency associated with this transition 
could, however, be misinterpreted as a radially decreasing accretion rate (dotted line). 



